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Abstract 

We find solutions of a gravity-Yang-Mills-Higgs theory in four dimensions that rep¬ 
resent asymptotic anti-de Sitter charged black holes with partial/full gauge symmetry 
breaking. We then apply the AdS/CFT correspondence to study the strong coupling 
regime of a 2 + 1 quantum field theory at temperature T and finite chemical poten¬ 
tial, which undergoes transitions to phases exhibiting the condensation of a composite 
charged vector operator below a critical temperature T c , presumably describing p+ip/p- 
wave superconductors. In the case of p + ip -wave superconductors the transitions are 
always of second order. But for p -wave superconductors we determine the existence 
of a critical value a c of the gravitational coupling (for fixed Higgs v.e.v. parameter 
rhw) beyond which the transitions become of first order. As a by-product, we show 
that the p -wave phase is energetically favored over the p + ip one, for any values of the 
parameters. We also find the ground state solutions corresponding to zero tempera¬ 
ture. Such states are described by domain wall geometries that interpolate between 
AdS^ spaces with different light velocities, and for a given rhw , they exist below a 
critical value of the coupling. The behavior of the order parameter as function of the 
gravitational coupling near the critical coupling suggests the presence of second order 
quantum phase transitions. We finally study the dependence of the solution on the 
Higgs coupling, and find the existence of a critical value beyond which no condensed 
solution is present. 


1 Introduction 

In recent years the application of AdS/CFT or more generally gauge/gravity correspondence 
CO PI 0 to the study of condensed matter physics has attracted a lot of attention, provid¬ 
ing in particular gravitational descriptions of systems exhibiting superconductor/superfluid 
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phases a 0- Since in condensed matter physics we are typically dealing with systems at 
finite charge density and temperature, in the context of the AdS/CFT correspondence the 
dual gravity descriptions should be given in terms of gravitational models with a negative 
cosmological constant which admit charged black holes as vacuum solutions. In fact, a 
charged black hole naturally introduces a charge density/chemical potential and tempera¬ 
ture in the quantum field theory (QFT) defined on the boundary using the gauge/gravity 
correspondence. This set-up allows in particular to study phase transitions and construct 
phase diagrams in parameter space. 

The simplest model is provided by an Einstcin-Maxwell theory coupled to a charged 
scalar field that, in the framework of the AdS/CFT correspondence, is dual to a scalar 
operator which carries the charge of a global U( 1) symmetry. It has been shown that a 
charged black hole solution, interpreted as the uncondensed phase, becomes unstable and 
develops scalar hair at low temperature breaking the 17(1) symmetry near the black hole 
horizon 0 0. This phenomenon in general may be interpreted as a second order phase 
transition between conductor and superconductor phases, interpretation that is supported 
by analyzing the behavior of the conductivity in these phases [5]. There were also studied 
vortex like solutions that describe type II holographic superconductors 0 0 US! and more 
recently, spatially anisotropic, abelian models of superconductors [Hi- 

Soon after these “s-wave” holographic superconductor models were introduced, holo¬ 
graphic superconductors models with vector hair, known as p -wave holographic supercon¬ 
ductors, were explored numerically first in [12] and [T3] (for a recent analytical treatment, 
see [33]). The simplest example of p-wave holographic superconductors may be provided 
by an Einstein-Yang-Mills theory with SU( 2) gauge group and no scalar fields, where the 
electromagnetic gauge symmetry is identified with an 17(1) subgroup of SU( 2). The other 
components of the SU( 2) gauge field play the role of charged fields dual to some vector 
operators whose non-zero expectation values break the 17(1) symmetry leading to a phase 
transition in the dual field theory. 

More recently, solutions to gravity-matter field equations where both scalar and vec¬ 
tor order parameters are present were considered; they describe systems where competi¬ 
tion/coexistence of different phases takes place [I5]-[T8]. 

Regular, self-gravitating dyonic solutions of the Einstein-Yang-Mills-Higgs (EYMH) equa¬ 
tions in the BPS limit and asymptotic to global AdS space, were constructed time ago in [0D3] 
and |2Uj . They were extended to dyonic black hole solutions in [21] and [22]; where they were 
interpreted as describing a so-called p + ip-wave superconductor (isotropic) system at finite 
temperature in the condensed phase. The purpose of the present work is to generalize pre¬ 
vious results by finding more general black hole solutions of EYMH in asymptotically AdS 4 
space with finite mass and electric charge density, to interpret them via the gauge/gravity 
duality as describing phases of a strongly coupled field theory, and to construct the corre¬ 
sponding phase diagrams 0. More specifically, in first term we start by revisiting the analysis 
of [21] [22] , verifying the existence of second order phase transitions all along the parameter 
space. It was found (see for example [23]) that some holographic systems pass from a second 
order phase transition as a function of the temperature in the non back-reaction limit to 

1 We will be considering the usual plane horizon ansatz, relevant to study condensed matter systems with 
translational invariance. Under these circumstances, the magnetic charge density of the dyon solutions in 
[21] [22] disappears. 


2 



a first order when the gravitational coupling exceeds a certain value. Such a phenomenon 
occurs in holographic superfluids when the velocity is high enough [24] [[25], and it was mea¬ 
sured in certain types of superconductors [26] [27] [28]. We have found this kind of behavior 
in our system in the anisotropic case, finding condensed solutions and constructing the phase 
diagram. Second, we compute free energies and find that for any set of values of the free 
parameters that determines the solutions, the anisotropic phase is energetically favored over 
the isotropic phase, as conjectured in other contexts [22] [3D] [31] . Third, we analyze the zero 
temperature limit, case that had not been addressed before; for low enough gravitational 
coupling we find solutions which spontaneously break the U( 1) symmetry and have zero 
entropy, and so describe the true ground state of the system. For gravitational couplings 
higher than a critical value the solution disappears, which is interpreted as a second order 
quantum phase transition. Lastly, we study the effect of a non zero Higgs potential on the 
system. 

The paper is organized as follows. In section 2 we present the model and write the 
translational invariant ansatz for the fields and the equations of motion that reduces to 
a a nonlinear system of coupled ordinary differential equations. In section 3 we present 
generalities of the systems to be studied at non-zero temperature, in particular the analysis 
of the holographic map to be used. In section 4 we present the numerical results concerning 
the “BPS limit”, i.e. null Higgs potential, including computations of free energies. Section 
5 is devoted to the study of the zero temperature case and the description of the ground 
state of the superconductor, including the presence of quantum phase transitions as function 
of the gravitational coupling and variable Higgs vacuum expectation value (vev) rh w . In 
section 6 the effect of a non-zero Higgs potential is considered. A summary and discussion 
of the results is given in section 7. Finally two appendices are added, one containing the 
boundary expansions of the fields and other containing the equations of motion and free 
energy in other parameterization commonly used in the literature. 

2 The gravity-Yang-Mills-Higgs system 

2.1 The model 

We consider a gravity-Yang-Mills-Higgs system in a 1 + 3 dimensional space-time with 
Minkowski signature (—b++). We take SU (2) as the gauge group, with generators satisfying 
the algebra, 

[Xai -W] = e a bc X c ; a, b, c — 0,1,2 , £012 = +1 (2-1) 

and the scalar field in the adjoint representation, H = H a X a . The full action to be consid¬ 
ered is, 

S = S (bulk) + S m) + S (ct) (2.2) 

where 

s(Mt) = L dl *^{hi R+ i?)~ih F -« F ° MN 

- i D M H a D M H a - j ( H a H a - Ho 2 ) 2 ] 
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S {GH) = 7 T 2 f d 3 x^\ 2 K (2.3) 

2« JdM 

where k, e and A are the gravitational, gauge and scalar couplings respectively, L is the AdS 
scale related to the negative cosmological constant through A = —3/L 2 , and H 0 > 0 defines 
the vacuum expectation value of the Higgs held (and so the boundary condition at infinity, 
see below in (I2.14p ). As it is well-known, the Gibbons-Hawking term S is necessary to 
have a well defined variational principle [32], where K = V a n a is the trace of the extrinsic 
curvature, and h and n the induced metric and normal vector on dXi. The counter-term 
action will be discussed in section QJ 

The held strength Fjfa 1 v and the covariant derivative Dm acting on the Higgs triplet H a 
are dehned as, 

F a MN = d M A a N - d N A a M + e a b c A b M A c n ; D M H a = d M H a + e a6c A b M H c (2.4) 

Let us consider coordinates (x^,y) and an ansatz preserving translational invariance in 
the coordinates {ok 1 , y, — 0,1, 2}, 

9 = ~ f(y) A(y ) 2 dx ° 2 + y 2 (c(y ) 2 dx 12 + dx 2 ^ + L 2 

A = L~ l ( dx° J(y) X 0 + dx 1 K x (y) A, + dx 2 K 2 (y ) X 2 ) 

H = H 0 H{y) X 0 (2.5) 


In what follows it will be convenient to introduce the dimensionless coupling constants, 


a = 


K 

e L 


rhw = eH 0 L 


An = e 2 HA L 4 A 


( 2 . 6 ) 


The gravity equations of motion (e.o.m.) derived from (I2.3[) result, 


(y f(v))' + 3 y 2 - y 2 f(y) 
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A'(y) A(y) (y 2 c'(y)Y 


A{y) 2 c(y) V A (y) 

{y 2 f{y)A{y)d{y))' 


A(y) c(y) 


while that the matter e.o.m. are, 
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\f(y)A(y) i 


m w 2 H(y)‘ 


(2.7) 
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l 2 


J\y)) = 


v, = 
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i MvfMvf 

2 2 y 2 c(y) 2 

(2.9) 


We will start by considering the “BPS limit” Ao = 0, but conserving the crucial Higgs 
vacuum value H 0 > 0. The effect of a finite Higgs coupling will be considered in section 6. 


2.2 Boundary conditions 

We will search for charged black hole solutions which present a horizon at y — where 
f(Vh) = 0. The associated Bekenstein-Hawking temperature of the black hole is given by, 


Tbh = TkL A( yh ) f'iVh) 

The ansatz (and e.o.rn.) are invariant under the scale transformations, 


(x°;A(y),J(y)) 
( a;1 ; c (y), Ki(y)) 


B‘ c(y), 8' Ki(y)^ 


( 2 . 10 ) 


( 2 . 11 ) 


They allow to fix some normalization imposing the b.c., A(y),c(y) v —m 1, in such a way 
that the x^’s are identified with the minkowskian coordinates of the boundary QFT, and 
(12.1 Op with its temperature. Furthermore there exists another scaling symmetry, 


M,y) —,7 y 


x h 
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f{y)^i 2 f{y) , Mv) i K i{y) , J{y)^iJ(y) 


( 2 . 12 ) 

that if yh ^ 0, allows to fix — 1 q Since now on we will fix the position of the horizon in 
this way, having in mind that we have to consider only scale invariants quantities. 

In [2T] and [22] solutions to (j2.7f) - fl2.8f) with a horizon and asymptotically AdS^ were 
studied. More specifically, there were found solutions with Ad = K 2 = K and the following 
boundary conditions; near the horizon y —y 1 + , 


f(y) = fi (y ~ !) + 0[(y - l) 2 ] 

2 This is the case except when we consider the zero temperature limit. When back-reaction is not taking 
into account yh = 0 corresponds to AdS space; when it is considered, yh = 0 is imposed in order to get a 
true ground state description and (12.1211 can be used to fix the chemical potential, see section [5] 


5 






















(2.13) 


A (y) 
c(y ) 
H(y) 
K(y) 

J(y) 

while on the boundary y —> oo, 


a o + a i {y ~ 1) + 0 [(y — 1)"] 
c 0 + ci (y - 1) + 0 [(y - l) 2 ] 
ho + hi (y — 1) + 0 [(y — 1)"] 
ko + h(y-l) + 0[(y - l) 2 ] 

ji (y -1) + C>[(y - l) 2 ] 


/(y) 

^(y) 

c(y) 

H(y) 
K(y ) 

■/(y) 


2 
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y 
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i + ... 

Hi 

1 H-y + ' 

yi 

y Kl 

, 0 + h + . 

y 


where consistency with the e.o.rn. and finiteness of /d(y) fixes to be, 


«i («i 


1) — 7Ti\y 2 



- +m w 2 


(2.14) 


(2.15) 


For more about the b.c. at the boundary, we refer the reader to the appendix A. We will 
adopt the b.c. (l2.13IWI2.14p in this paper except in section [5] where the b.c. on the horizon 
will have to be modified. 

The bulk theory is invariant under the gauge group S77(2); however the b.c. on the Higgs 
field, H{y) ? 1 , breaks this invariance to the U{ 1) generated by X 0 . With respect to 
this gauge subgroup the electric charge density of a solution is defined as usual by, 



J_ c(y) 
L 2 A{y ) 


y“ J (y) ly^-oo 


A 

L 2 


(2.16) 


As we show in section [4] at fixed couplings (a,rhw) a general solution to (I2.7p ~ fl2.8p with 
the b.c. (I2.13p ~ fl2.14p is determined by Jo, which is related to the U(l) chemical potential 
by, 

II = A“(oo) = J (2.17) 

From (I2.16p and (I2.17P the standard asymptotic expansion follows, 


A °o(y) 



y 


(2.18) 


Along this paper we will adopt /i as our scale. From (12.121) the dimensionless, scale invariant 
temperature is, 


Tbh _ o-o fi 
jJL 4 7T Jo 


(2.19) 
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where a 0 and /i are defined in (j2.13j) . A solution is determined by the three free parameters 
(a,rhw, J 0 ), and so the temperature (through the coefficients ao, /i) results a function of 
them □. 

In the analytic solution to the equations (I2.7p - fl2.8p that preserves the U(l)x 0 symmetry 
matter fields take the form, 

J(y) = Jo + - = Jo(l-~) ; Ki(y) = 0 ; H(y) = 1 (2.20) 

y V vJ 

where we imposed smooth behavior of the gauge field at the horizon which yields the con¬ 
dition J( 1) = 0, see the last line in (I2.13p . and then fixed J\ = — Jo 0. In what the metric 
functions concern, they correspond to the AdS Reissner-Nordstrdm (AdS-RN) black hole, 


m 

f(y) 


with temperature, 


c(y) = 



y -1 


y 2 


l 


i + 




1 a 2 Jo 2 1 
y + ^2~^ 

+ (y 2 + 2y + 3) (y - 1)) 


r - _J_ A _ AZl 

4ttJ 0 V 2 

The extremal, zero temperature AdS-RN black hole is defined by the relation a 2 Jo 2 


( 2 . 21 ) 


( 2 . 22 ) 
= 6 . 


3 Solutions at T > 0: superconducting state 

When the “magnetic part” of the gauge field is non-trivial, i.e. Ki(y) ^ 0 for some i = 1,2, 
the solution breaks not only the U(l)x 0 invariance, but also the invariance under rotations 
in the (x 1 , x 2 )-plane. According to the AdS/CFT dictionary this hair is interpreted as 
a spontaneous breaking of a global 17(1) symmetry present in the boundary QFT, whose 
currents take an expectation value, 

{J?{x))~KiS? ; *> = 1,2 (3.1) 

Giving that the order parameter is dual to (a component of) the gauge field we are pre¬ 
sumably modeling a p- wave superconductor |T3] . The normal state of the superconductor is 
described by the AdS-RN solution (I2.20p - (l2.2ip ; such solution is energetically favored until 
a critical temperature T c is reached; when T < T c the non symmetric, hairy solution gives 
rise to a superconductor phase. 

We remark that with the b.c. on the Higgs field we are breaking explicitly the gauge 
group from SU{ 2 ) to U(l)x 0 ] this yields a mass for the “W” gauge bosons, 

m w = e H 0 (3.2) 

3 In EYM systems where the Higgs field is not present the temperature is function of the only free 
parameter of the theory, a, see [23] . 

4 When yh = 0, J(y) = Jo is just the chemical potential and the metric solution is AdS space; it describes 
the uncondensed phase when the temperature is zero, see section 5. 
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The problem is thus the following: can we find under this condition a solution with K,(y) ^ 0 
that breaks spontaneously the U( l)x 0 ? In the boundary QFT this is then interpreted as the 
breaking of a global {7(1) symmetry as it happens in superfluids and superconductors with 
weakly coupled photons. From here we identify T c with the critical temperature of the phase 
transition in the QFT. 

We will consider two cases. 

• The isotropic case: K(y) = K\ (y) = K 2 (y) 

Although both gauge and rotational symmetries are broken by a hairy solution, a 
configuration (12.5p with K\ (y) = K 2 (y ) preserves the diagonal subgroup, 

({/(l) A -„ x SO(2U) Jiag , fact that is manifest in (18.ip [12] , This configuration give 
rise to an energy-momentum tensor isotropic in the x 1 -x 2 plane; therefore the metric 
function c(y) must be a constant, even when back-reaction is taken into account. 

This kind of configurations were first studied in |21j-[22j. using relaxation methods. 
We will re-obtain these solutions here for later use by using shooting methods. 

• The anisotropic case: K(y ) = K\(y) ; K 2 (y) = 0 

As stated above, a configuration with K\ = 0 preserves the U(l)x 0 and spatial rota¬ 
tions. When K\ (y) develops a non zero value the gauge symmetry U(l)x 0 breaks, and 
the condensate K\ (y) X\ dx 1 choose a direction x 1 as a special one. Then if we take 
into account back-reaction effects the system cannot support the condition gn = g 22 
[13] , Due to this fact T x i x i ^ T X 2 X 2 and the function c(y) can not be a constant; in 
conclusion the system will be in an anisotropic phase. 

In both cases the vacuum expectation value in the d = 3 field theory of the current 
operator Ok, dual to the function K associated with the magnetic field in the bulk, follows 
from the identification, (Ok) ~ K\ with K\ defined in (12.14h : K\ = K\ (T) can be taken 
as the order parameter that describes the phase transition of the system. As discussed for 
different models 0-|22| one can interpret this result by stating that a condensate is formed 
above a black hole horizon because of a balance of gravitational and electrostatic forces. 
From the asymptotic behavior in (12.Iff we get the dimension A [Ok\ of the operator Ok 

ra 

AlOjf] - 1 - x - ) - (v'l - Ani v r (3.3) 

From numerical solutions we conclude that a finite temperature continuous symmetry 
breaking transition takes place so that the system condenses at a critical temperature T c , as 
can be seen from the behavior of K\ (T) for T « T c in figures 3, 4 and 5. Furthermore, we 
compare the free energies corresponding to both phases in figures 8 and 9, finding that the 
anisotropic phase is favored, see f29j [30] [31] for related results. 

4 Numerical Solutions 

We analyzed numerically equations (12.7[) ~ (j2.8p and found solutions that satisfy the required 
b.c. (I2.13p - (l2.14p in a wide region of the parameter space, that lead to the phase diagram 
in figure 1. Such solutions in the anisotropic case are shown in figure 2. 


















Before presenting the results, we think is worth to spend a few words on the method used. 
As discussed in appendix A, after fixing some normalization and asking for finiteness the 
solution near the boundary admits the expansions in equation (1A.1I) . and is determined by six 
constants, (Fi, C±, Jo, Ji, Ad, Hi). However the b.c. on the horizon impose five conditions. 
The first two come from the definition of the horizon and the regularity of the gauge field, 


/Cl) = 0 ; 7(1) = 0 


(4.1) 


They essentially fix the mass (~ Fi) and the charge density (~ Ji) of the black hole. The 
remaining three conditions fix (Ci, K\, Hi) and are obtained from an analysis of the (singular) 
behavior of the e.o.m. near the horizon, 


c'(l) 
^(1) 
H\ 1 ) 


2 2 C (l) 

a m w 


Ki (1) 

m 

H( 1 ) 


KiOl 

m v c(i)2 


+ A" 2 (l) 2 ) H( 1) 


(K 2 (l) 2 + m w 2 H(l) 2 ) 


/'(l) V c(l) 


Ad(l ) 2 


+ k 2 (i)‘ 


(4.2) 


Therefore the only additional free parameter that determines the solution is Jo, i.e. the 
chemical potential (12.171) . In practice we integrate the system from the horizon, where 
according to (I4.1l) - (l4.2p the free parameters are, 


</'(!)= ji ; K(l) = ko ; H(l) = h 0 ; A(l) = a 0 ; c(l) = c 0 (4.3) 


as dehned in (12.13ft . These parameters are selected in such a way that the solution matches 
the conditions on the boundary (12.14[) . 


A(oo) = c(oo) = H(oo ) = 1 ; K( oo) = 0 ; J(oo) = J 0 (4.4) 



Figure 1: Phase diagrams in the isotropic case (left) for rhw = 0.1 (blue), rhw = 0.4 (green), 
and in the anisotropic case, for rhw = 0.1 (center) and rhw = 0.4 (right). 


Figure 1 displays the phase diagrams in the a — T plane, for two different values of rhw 
In the white regions only the normal or uncondensed phase is present. In the isotropic case 
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the system experiments second order phase transitions along the blue (rhw = 0.1) and green 
(■rhw = 0.4) curves. In the anisotropic case, in the blue and red regions the condensed phase 
is the thermodynamically preferred phase. The blue line until the black point indicates a 
critical line of second order transitions. The black point signals the coupling cq c beyond 
which the transitions become of first order along the red line, while that the blue and green 
lines that continue after this tri-critical point represent spinodal lines. The critical red curve 
of first order phase transitions ends in the red point at T = 0, which represents a quantum 
phase transition, signaling a critical coupling a above which the condensed phase ceases to 
exist, see Section [5j By comparing both graphs we can see that both a tc and a decrease 
with increasing fri\y. Similar phase diagrams were obtained in reference [34] in absence of 
Higgs fields. 

In figure 2 the fields are shown as functions of the coordinate y, at fixed J 0 and rhw and 
for different a’s. For a c ~ 0.8825 a second horizon appears, as displayed from the curves 
corresponding to f(y)/y 2 . The uncondensed and condensed phases are separated by a curve 
on which the formation of the second horizon takes place for a given critical temperature 
determined by the gauge boson mass rhw and J 0 . In the isotropic case the curve is displayed 
in figure 1 for two different values of rhw (blues and green lines) and it coincides with the 
critical curve on which the phase transitions take place. On the other hand, in the anisotropic 
case the curve coincides with the critical curve (blue line in figure 1) until the tri-critical 
point atc, and it continues through the spinodal curve in green. 

In figure 3 the order parameter Ki in the isotropic case is plotted as a function of the 
temperature for different values of rh w , at fixed gravitational coupling a = 0.7. In this case 
the transition is of second order independently of a, in agreement with fT2j . 

Figures 4 and 5 shows the order parameter K\ as function of T in the anisotropic case 
from two perspectives: at fixed rhw — 0.4 and varying a in figure 4 and at fixed a = 0.7 
and varying rhw hi figure 5. From figure 4 it is seen that for a tc ~ 0.53 K\ becomes 
multi-valued, fact that signals the passage from second to first order phase transitions as 
corroborated from the free energy computations of the next subsection. This phenomenon 
has been found recently in p -wave superfluids by studying the role of the back-reaction 
in the phase transitions [23] (for experimental results on first order phase transitions in 
superconductors, see [26] [2?] [128]). By comparing figures 3 and 5 it is observed that the 
temperature at which the order parameter becomes zero is the same in both cases, and 
that the critical temperature decreases when rh w increases, what can be interpreted as the 
presence of the Higgs held hinders the condensation. Furthermore, we have checked that 
near T c and for weak gravitational couplings a < atc, K\ behaves like (T c —T) 2 , indicating a 
second order phase transition with mean held exponent | as usually happens in holographic 
descriptions of critical systems in the limit of large number degrees of freedom. 
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f(y)/y 2 
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Figure 2: Solutions for the fields (f(y),A(y),c(y),K(y),J(y),H(y)) with the b.c. (12.13H - 
(12.14ft in the anisotropic case. The curves correspond to Jo = 6, rhw = 0.4 with different 
fixed values of a = 0.0 (black), 0.1 (green), 0.2 (light red), 0.3 (light blue), 0.4 (fuchsia), 0.5 
(brown), 0.6 (gray), 0.7 (orange), 0.8 (red), 0.88 (blue). One can appreciate from the curves 
corresponding to f(y)/y 2 the formation of a second horizon at a c ~ 0.8825. The analogous 
solutions for the isotropic case can be found in [22] . 
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Figure 3: The order parameter Ki/J 0 1+K1 = (Ox 1 )/Jo 1+K1 is plotted in the isotropic 
case at fixed a = 0.7, for different values of rh w = 0.01 (black), 0.1 (green), 0.2 
(orange), 0.3 (brown), 0.4 (blue) that correspond to the critical temperatures T c = 
0.017056,0.016371,0.014174,0.011215,0.00791 respectively. Near the critical temperature 
the order parameter behaves like (T c — T) 1 / 2 . 


AW 1+1 



Figure 4: The order parameter K x /J q 1+K1 = {Ok x )/Jq 1+K1 is plotted in the anisotropic case 
at fixed rhw — 0.4, for different values of a = 0.50 (black), 0.53 (blue), 0.55 (red). In the 
inset is displayed the multi-valuation of the order parameter for a tc ~ 0.53. 
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Figure 5: The order parameter Ki/J 0 1+K1 = {Ok x ) / Jo 1+K1 is plotted in the anisotropic case 
at fixed a = 0.7, showing the phase transition at different values of rh\y = 0.01 (black), 0.1 
(green), 0.2 (orange), 0.3 (brown), 0.4 (blue), 0.32 (blue). 


4.1 The free energy 

According to the AdS/CFT correspondence, the free energy of the QFT is given by, 


F = T Spiirl = 


dy I dx L euc i 


(4.5) 


'Vh 


From (12.3ft and using the e.o.m. the bulk contribution to the free energy density can be 
written as, 

i rv°° 


f {b 


mlk) _ 


2 e 2 L 3 


dyA(y) c(y)y 
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j'{yy 


c(yY 


6 Ao 


a 2 2 
K^yfK^yf 


(H(y) 2 -iY 
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A(y) 2 y 2 f(y)A(y ) 2 \ c(y ) 2 

The Gibbons-Hawking contribution is, 


f {GH) = 


2 e 2 L 3 


a- 


fiy)* (y 2 f(y)~ 2 My)c(y) 


(4.6) 


(4.7) 


Here we have introduced yoo to regularize the expressions since they present divergent terms. 
To this end we introduce a counter-term action 


g{ct) — 


1 


2 K 2 


[ d 3 x y/\h\ —- 
IdM L 


which give rise to the following contribution to the free energy density, 


/ <ct) = 


4 


2 e 2 L 3 a 2 


— (f{y) 1 / 2 A{y)c(y)y 2 ) 


Uoo 


(4.8) 


(4.9) 
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The total free energy density of the system / is then given by, 

/= lim (/( Mfe > + / (aff ) + / (ct) ) (4.10) 

J/oo—>00 V 7 

We remark that in order to analyze the results the right thing to do is to work with the scale 
invariant free energy density, 

Figures 6 and 7 show the evolution of the free energy density (14.111) with the mass of the 
gauge boson for two different values of a, in the isotropic and anisotropic cases respectively. 
Figure 6 displays the continuity of / at the critical temperature (where the free energy density 
of the uncondensed phase intersects the curve of the condensed phase) for both values of a, 
for any rhw, fact that indicates the second order character of the phase transition as the 
behavior of K i in figure 3 suggested. In figure 7 instead it is observed the discontinuity in 
the first derivative of the free energy density at the critical temperature for a = 0.7 > a tc 
for any rhw, signaling a first order phase transition. In both cases the critical temperature 
decreases with growing rhw, in agreement with the analysis of the behavior of the order 
parameter made above. 

In figures 8 and 9 the free energy densities of the isotropic and anisotropic phases are 
compared for two values of the gravitational coupling, a < a tc (figure 8) and a > a tc (figure 
9). From them one can see that the free energy density of the anisotropic phase, no matter 
the region where the value of a is, i.e. if first or second order phase transitions take place, 
is lower than the free energy density of the isotropic phase. That is, the anisotropic phase 
is always energetically favored over the isotropic one. 




Figure 6: The free energy density / is plotted as function of the temperature in the isotropic 
case for two values of the gravitational coupling constant, a = 0.4 (left) and a = 0.7 (right), 
at different values of rhw = 0.01 (black), 0.1 (green), 0.2 (orange), 0.3 (brown), 0.4 (blue). 
The red curve represents the free energy density of the uncondensed phase. 
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Figure 7: The free energy density / is plotted as function of the temperature in the 
anisotropic case for two values of the gravitational coupling constant, a = 0.4 (left) and 
a = 0.7 (right), at different values of rh w = 0.01 (black), 0.1 (green), 0.2 (orange), 0.3 
(brown), 0.4 (blue). The red curve represents the free energy density of the uncondensed 
phase. 
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- 0.11 
- 0.12 
-0.13 
-0.14 
T -0.15 



Figure 8: The free energy densities for the isotropic (black) and anisotropic (blue) cases 
are plotted as function of the temperature for different values of rhw, at fixed a = 0.4. The 
red curve represents the free energy density of the uncondensed phase. 
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Figure 9: The free energy densities for the isotropic (black) and anisotropic (blue) cases 
are plotted as function of the temperature for different values of m w , at fixed a = 0.7. The 
red curve represents the free energy density of the uncondensed phase. 


5 Zero temperature solutions 


In this section we will address the problem of quantum phase transitions in three dimensional 
p -wave superconductors (anisotropic case), i.e. transitions at T = 0, that to our knowledge 
was not considered before in the literature (see however [37] [38] [39] [40] for related studies 
in other settings). 

It is known that when a charged AdS black hole is driven to a state of zero temperature 
it becomes extremal, but its entropy is different from zero and then it can not describe the 
ground state of the superconductor that we are presumably modeling holographically. To 
reach our goal the radius of the black hole needs to become null, to agree with the third 
law of thermodynamics and really describe the quantum ground state [29] [4Jjj. So we must 
impose that yh = 0, i.e. the coordinate y E [0, oo). A very important thing from a technical 
point of view is that while the asymptotic behavior of the fields is as in (I2.14jh the expansions 
near the horizon drastically change with respect to the T > 0 case. At leading order the 
(non analytical) behavior of the fields for y —> 0 + is, 


f(y) 


~ ^ _ 2fc 0 

2 2 k ojo 2 e~ £ o« 

y — a --- 

2 c 0 a 0 2 y 
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fcp 
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(5.1) 


The independent constants are k 0 , jo, ho,do,co- Such constants are chosen in the same way as 
in the T > 0 case, see (I I. Hi . From (15. ip it follows that near the horizon the solution is another 
AdS 4 space. The solutions that describe the quantum ground state of the superconductor 
in the condensed phase are therefore domain walls interpolating AdS 4 spaces with the same 
radius L but different light velocities in both directions, in virtue of the fact that do ^ 1 and 
Co 7 ^ 1. More explicitly, 


uv 

Vi _ Vh 

c{yh) 

/ / (2/00) 

Mv 00) 

_ c 0 

v RV _ yh 

/ fiVoo) 

A{yoo) _ 1 

v[ R 2/00 

c(y oc) 

V f(yh) 

MVh) 

do 

V{ R 2 /oo \ 

/ f{yh) 

MVh) do 


(5.2) 


On the other hand, the uncondensed phase is described strictly by AdS 4 space and J(y) = J 0 , 
which replace the AdS-RN solution (I2.20p - (l2.2ip . Interestingly, we found that above a certain 
a the solution representing the condensed phase disappears and the only solution that exists 
is AdS space. This result can be guessed from the following analysis borrowed from [|38] 
(see also [23]). At very low temperatures the normal phase is nearly represented by the 
extremal, zero temperature Reissner-Nordstrom solution (12.20112. 22]) with Jo 2 = Jj, whose 
near horizon geometry is AdS 2 x 9ft 2 . If we perturb this solution with a non-zero gauge field 
K\ (y) = K(y ), from the first equation in (12. 8p its linear equation in this background results, 


0 = (d 2 d 2 p + 2 pd p - rh eff 2 ) K(y) 

where p = y — 1, which is just the wave equation for AdS^ with an effective mass, 


^e// — 


1 

6 



(5.3) 


(5.4) 


So, the instability to form SU (2) vector hair at low temperature is just that of scalar fields 
below the BF bound for AdS 2 , m 2 BF = —That is, when m e ff 2 < rh 2 BF one could wait 
that AdS vacuum gets unstable and the system prefers to be in the phase described by the 
superconducting black hole solution with non abelian hair. Thus we get a plausible condition 
for instability, 


ry <C (y — - 

\ ^q uess — 3 , o 

y f + m w 2 


(5.5) 


Figure 10 shows the fields for different values of the gravitational coupling 0. In the 
example showed (rhw = 0.4) we obtain a ~ 0.8825 ~ cx c |r> 0 ; see figures 2 and 10. 


5 We use the scaling symmetry (12.121) to fix J 0 = 1. 
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Figure 10: We display for T = 0 the solutions f(y)/y 2 , A(y), c(y ), H(y), K(y ) and J(y) as 
functions of the coordinate y from the horizon, now located at y — 0, towards the boundary, 
for rhw = 0.4 and different values of a = 0.0 (blue), 0.3 (orange), 0.6 (green), 0.8 (fuchsia). 
One can appreciate from the curves corresponding to f(y)/y 2 the formation of a second 
horizon when a approaches a ~ ck c |t>o ~ 0.8825; for larger values of a the asymmetric 
solution ceases to exist. 


We worked out the solutions for different values of the parameter rhw, in Table I the 
corresponding critical couplings are shown. It is observed that a decreases with growing 
and that a guess < <5, what is consistent with the instability analysis made before. 

On the other hand, in Figure 11 it is shown the order parameter as a function of the 
coupling a. We have verified that the behavior near the critical coupling is of the type, 

A0(a) ~ (b — a)? (5.6) 

consistent with the existence of a second order phase transition in the mean field, large 
number of degrees of freedom limit. 
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rh w 

a 

0.1 

1.03 

0.2 

0.97 

0.3 

0.94 

0.4 

0.88 

0.5 

0.86 

0.6 

0.84 

0.7 

0.80 

0.8 

0.76 

0.9 

0.74 

1.0 

0.70 


Table 1: Critical gravitational couplings a for different values of rhw at T = 0. 
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Figure 11: The order parameter K\ is plotted as function of a at T = 0, for different values 
of rh w — 0.5 (blue), 0.6 (lightblue), 0.7 (green), 0.8 (orange), 0.9 (brown), 1.0 (red). The 
values for which = 0 define the critical couplings a. 


6 Analysis for A^0 


In this section we will study the effect of a non-zero Higgs potential as introduced in (12.31) . 
specified by the Higgs vev scale Hq and the strength A. For simplicity we will work in the 
no back-reaction limit a = 0, although the new insights does not depend on this fact. 

In the conventions of Appendix B, the e.o.rn. (j2.8j) in the anisotropic case reduce to, 
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(. f{n)K\u ))' 
2 I f ( U ) 
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H'(u)\ = 


= m w 


H[u) 2 J(uf 




/(«) 


K(u) 


w n2 Ao H(u) 2 — 1 . TT . . 
K(u) 2 + ^ - H{u) 


mw 
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J"{u) = 


K(u) 

f(u) 


J{u) 


( 6 . 1 ) 
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where f(u) — 1 — u 3 . 

The existence of a non-zero A does not modify the behavior of K(y) and J(y) on the 
boundary that remain as in (12.14[) . but it does in the Higgs case where we now have, 

H(u) = 1 + H^u a - +--- + H+ u A + + ... (6.2) 

where, for general Aq and mw, 


A ± = 



(6.3) 


A hrst well-known fact is that reality of A± necessarily implies the BF bound Ao > — | rhw 2 - 
When we are in the window — § < < — f both modes are normalizable and lead to 

consistent quantization and we can impose = 0 or H + = 0. If Ao > — § mw 2 , the 
condition //_ = 0 must be imposed 021 . We will consider for definiteness the case Ao > 0. 

A very interesting fact is that, besides the existence of a bound from below for the Higgs 
coupling as stated above, a straight analysis of the solution near the boundary u = 0 yields 
the result that a bound from above is also present. We find that there exists a critical value 
Aq defined by, 

-^3=2 + (k 1 -1) (2/ci + 3) (6.4) 

mw 

such that for A 0 > K the condensed solution ceases to exist. In the example considered 
below rhw 2 = 1, Aq ~ 5 .854. This is so for both the isotropic and anisotropic cases. 


w“ w * 1 



Figure 12: The condensate K\ as a function of the temperature T, for mw 2 = 1 and Ao = 0 
(black), 0.25 (orange), 0.5 (brown),0.75 (blue), in the anisotropic (left) and isotropic (right) 
cases. 
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Figure 13: The free energy / as a function of the temperature T for m\y 2 = 1 and Ao = 0 
(black), 0.25 (orange), 0.5 (brown), 0.75 (blue) in the anisotropic (left) and isotropic (right) 
cases. 


/ / / / 



Figure 14: Comparison between the free energies of the anisotropic (blue) and isotropic 
(black) cases for mw 2 = 1 and, from left to right, A 0 = 0, 0.25, 0.5, 0.75. 
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Figure 15: The free energy as a function of the temperature in the anisotropic (left) and 
isotropic (right) cases, for fixed A 0 = 0.25 and different values of rhw — 0.5 (blue), 0.7 
(brown) and 1.0 (black). The red curve represents the free energy of the normal phase. 


Figures 12 and 13 show the condensate and the free energy respectively as functions of the 
temperature, for a fixed rhw and different Higgs couplings l|. One can appreciate that both 
the order parameter and the free energy decreases with increasing strength of the potential 
Ao; however the critical temperature does not change and the phase transitions remain of 
second order. 

In Figure 14 the free energies of both isotropic and anisotropic cases as functions of the 
temperature are plotted together for comparison, for fixed Higgs scale and different Higgs 
couplings. The anisotropic phase always remains energetically favored over the isotropic one. 

As stated before, a critical value (16. 4 j) above which the condensed phase does not exist is 
present, and it was verified by our numerical calculations; such result occurs in both isotropic 
and anisotropic cases. 

Figures 15 show the free energies in function of the temperature for a fixed Ao, for different 
gauge boson masses. It is observed that they increase with increasing rhw■ The curves are 
similar to the left curves in figures 6 and 7. 

7 Conclusions and outlook 

In this paper we have investigated four dimensional solutions of black holes with non-abelian, 
SU (2) hair introduced by Yang-Mills gauge bosons and a non trivial Higgs field in the adjoint 
representation, whose v.e.v. triggers the breaking of the gauge symmetry to a 77(1) subgroup 
under which the black hole is charged. 

In the spirit of the AdS/CFT correspondence, the symmetric solution given by the AdS- 
RN black hole when the temperature is positive and AdS when T = 0, describes the uncon¬ 
densed phase of the dual three dimensional QFT. A solution with non-abelian hair generically 
breaks the 7/(1) gauge symmetry together with the rotational symmetry, and is interpreted 
as describing a condensed phase of the QFT. The order parameter is the coefficient of the 

6 For higher Higgs strengths towards the critical value the curves does not experiment significative changes. 
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leading order term of the magnetic component of the gauge field, and thus the systems de¬ 
scribed are generically termed p -wave superfluids/superconductors. We have considered two 
cases. The isotropic case that describes p + ip-wave superconductors where the diagonal 
subgroup of U[l)gauge x SO(2) rot is preserved, and the anisotropic case where no symmetry 
is preserved. In both cases we get phase transitions at critical temperatures that decrease 
when the gravitational coupling grows and in the case of anisotropic superconductors the 
phase transitions become of first order for large gravitational couplings [23]. These results 
are summarized in the phase diagrams presented in figure 1. 

We also find solutions that describe the zero entropy ground state of the p -wave super¬ 
conductor, showing the existence of phase transitions from the normal phase (described by 
AdS space) to this condensed phase, that is present below a certain value of the gravita¬ 
tional coupling a. These transitions are of second order, according to the behavior (15.6|) of 
the order parameter near the transition obtained from figure 11. Such states are described 
by domain wall geometries that interpolate two AdS spaces. The occurrence of AdS space 
near the horizon, with the the same scale as the AdS in the boundary but different light 
velocities, presumably indicates that there is an emergent scale invariance in the T = 0 limit 

[23] [21]- 

Finally we study the effect of considering a non zero Higgs potential, ft was found that for 
Higgs coupling constants greater than a critical value A c , the solution collapses to the normal 
one. This fact relies on the ultraviolet behavior of the system; in particular is independent 
if the back-reaction is considered or not. For A below A c the system and its thermodynamic 
variables behave qualitatively as in the case A = 0, but with lower free energy. 

A very relevant fact conjectured in the literature [12] [13] [293 f° r systems without Higgs 
fields that we have explicitly addressed in this paper including them together with the 
corresponding Higgs potential, is that below the critical temperature in all the parameter 
space we found that the free energy density of the anisotropic solution is lower than that of 
the isotropic one, indicating that the p -wave superconductor phase is more stable that the 
one corresponding to the p + ip-wave superconductor. This result is illustrated in figures 8, 
9 and 14. 

We believe it is worth to make the following remarks. It is straight to see that if we 
switch off the Higgs field the e.o.m. of the p -wave superconductor are recovered. However if 
we switch off the magnetic part of the gauge field, Ki = K 2 = 0, it is seen that we do not 
recover the e.o.m. of a s-wave superconductor. This is due to the fact that we are switching 
the Higgs field in the Ao-direction, not in the Xi-direction. This lead us to conclude that 
in our set-up the Higgs field can not condense spontaneously since the temporal component 
of the gauge field (which plays a fundamental role in the condensation) is null. Therefore 
we will never have competition between s-wave and p-wave phases, as it takes place in the 
cases analyzed in references ra-ra where the matter field ansatz is slightly different and 
the vev Higgs field is put to zero. Furthermore it is not difficult to see from the e.o.m. that a 
configuration where a vev (13.11) is present necessary implies a non-trivial Higgs field; however 
the vev of its dual scalar operator 0 (x ), 

(0(x)) ~ H, (7.1) 

does not indicate any spontaneous breaking in view of the presence of the source H 0 . 
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We stress that, although they share some similar characteristics, the presence of the 
Higgs fields with the non trivial b.c. |/J(oo)| = Hq > 0, introduces a scale that makes our 
systems different from those considered in precedence since pT3j |37j . in which the scalars 
were not present. On one hand from the obvious fact that the system is larger and more 
complex; in particular we have three free parameters (a,mw, Jo) an d, among other things, 
the dimension (13. 3 p of the order parameter remains arbitrary. Instead, in EYM systems 
where the Higgs field is not present, the temperature for example is a function of just one 
parameter a = ^ [23]. On the other hand and more important, from the QFT point of view 
to which the systems we have considered are presumed to be holographically dual. This fact 
can be elucidate by studying the transport properties of the system, i.e. the conductivities. 
Even ignoring back-reaction effects, we finish with a system of fifteen coupled second order 
equations that results much more cumbersome to disentangle than in the Abelian case or in 
the absence of Higgs fields. We hope to report results in this direction in a near future [43] • 
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A Boundary expansions 


Along this paper we used shooting methods to get the solutions to the e.o.m. We present 
here the (next to) leading order behavior of the fields near the boundary that is necessary 
to carry out the numerics. 

For large y —> oo the fields admit the expansions, 
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where A,(0) = 1, a* > 0, etc., for i — 1, 2,.... The constants (Fi, Ci, J 0 , A, Ki, Hi) are free, 
all the other ones as well the powers (including those make explicit in (1A.1I) ) are determined 


24 







by the e.o.m. Q In the isotropic case, K(y) = Ki(y) = K 2 {y ), C(y) = 1 (C\ = 0) , they are 
given by, 
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0 < rh w 2 < 2 
mp/ 2 = 2 


A J Jo 


(2 Kl + l) (l + Kl) J\ 
K 1 2 Jo , Ai 
21 Jj ' 7 

_hi 

(2 Kl —1) (1 + Kl) Al 

All a 
9 Hi 2 


-j- < rhw 2 < 00 
0 < rhw 2 < x 

™ 2 _ 15 

— — 

2 < rhw 2 < 00 
0 < rhw 2 < 2 

mp/ 2 = 2 


(A.2) 


while that in the anisotropic case, K(y) = Ki(y), K 2 (y ) = 0, they result, 
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(A.4) 


B AH conventions 

In this appendix we write the e.o.m. in the conventions of references ESI. H- Besides being 
often present in literature, they proved to be convenient in some numerical computations. 
The ansatz is written as, 


9 = 


T 2 / Xr2 ] 2 

^ s / >2 7.2 , , / \2 7-2 , au 

+ g(u) dy + 
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A = dt J(u ) X 0 + da; K\ {u) X\ + dy K 2 (u ) X 2 
H = H 0 H(u)X 0 

The relation with the conventions used in the bulk of the paper are, 
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‘ The boundary conditions at the horizon leave just one free parameter (that we take J 0 ) that determines 
completely the solution, see Section [I] 
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J(y) = J(u) , K 1 (y) = K 1 (u) , K 2 (y) = K 2 (u ) , H{y) = H{u) (B.2) 

The gravity e.o.m. result (out the tildes), 
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where (1), (2), (3) are the r.h.s.’s of (12.71) written in the variables (IB.2ft . while that the matter 
e.o.m. are, 
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Finally, the contributions to the free energy density are, 
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